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Following our previous work [J. Ieda, T. Miyakawa, M. Wadati, cond-mat/0404569] on a novel 
integrable model describing soliton dynamics of an F = 1 spinor Bose-Einstein condensate, we 
discuss in detail the properties of the multi-component system with spin-exchange interactions. 
The exact multiple bright soliton solutions are obtained for the system where the mean-field 
interaction is attractive (co < 0) and the spin-exchange interaction is ferromagnetic (c2 < 0). A 
complete classification of the one-soliton solution with respect to the spin states and an explicit 
formula of the two-soliton solution are presented. For solitons in polar state, there exists a 
variety of different shaped solutions including twin peaks. We show that a "singlet pair" density 
can be used to distinguish those energetically degenerate solitons. We also analyze collisional 
effects between solitons in the same or different spin state(s) by computing the asymptotic 
forms of their initial and final states. The result reveals that it is possible to manipulate the 
spin dynamics by controlling the parameters of colliding solitons. 
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1. Introduction 

Solitons are a universal feature of nonlinear phenom- 
ena 1 in many different physical areas such as particle 
physics, nonlinear optics, plasma physics, fluid dynamics, 
and condensed matter physics. Recently, two teams in 
the US 2 and France 3 performed experiments that demon- 
strated matter-wave bright solitons with Bose-Einstein 
condensates (BECs) of gaseous 7 Li atoms. In general, 
solitons are formed under the balance between nonlinear- 
ity and dispersion. For atomic BECs whose macroscopic 
wave functions obey the Gross-Pitaevskii (GP) equation 
or nonlinear Schrodinger (NLS) equation, the former is 
attributed to the interatomic interactions while the latter 
comes from the kinetic energy. 

Either dark or bright solitons are allowable depending 
on the positive or negative sign of the interatomic cou- 
pling constants, respectively. Bright solitons created in 
the experiments are themselves condensates and prop- 
agated over much larger distances than dark solitons 
which, on the other hand, can only exist as notches 
or holes within the condensate itself. 4 Those observed 
matter- waves behaved like scalar fields since the spins of 
condensed atoms are frozen under additional magnetic 
fields. In line with this, many theoretical studies on the 
bright soliton formation and propagation of attractive 
BECs have been done mainly in the single component 
systems. 5-8 

Multi-component generalization of the soliton dynam- 
ics is very natural in the context of atomic BECs because 
of several ways to create such systems, i.e., the mixtures 
with two different atomic species/hyperfine states 9 and 
the internal degrees of freedom liberated under an optical 
trap. 10 Moreover, when one explores future applications 
of matter- wave solitons in atom optics, 11 such as atom 
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laser, atom interferometry, and coherent atom transport, 
the multiple condensates have various potential utilities, 
e. g., the multi-channel signals and their switching. 

So far, multiple solitons were theoretically investigated 
in two-component condensates 12 and also in an arbitrary 
number components system; 13 the nonlinearity in both 
cases consists only of the intensities, i.e., the squared 
absolute values of the components. The coupled NLS 
equations with such nonlinearity, often referred to as the 
intensity-coupled type, are extensively used for model- 
ing multi-mode optical solitons. 14 In contrast, the so- 
called spinor condensates, 15-20 being another candidate 
that supports multiple matter-wave solitons, have non- 
trivial nonlinear terms reflecting the SU (2) symmetry of 
the spins. The spin-exchange interactions that are the 
sources of the spin-mixing within condensates 21 ' 22 are 
the exceptions of the intensity-coupled nonlinearity and 
there is no analogue in conventional optics. 

In our previous work, 23 we discovered a novel inte- 
grable model which describes the dynamics of a spinor 
BEC and derived the multiple TV-soliton solution. We 
assume specifically an atomic condensate in the F = 1 
hyperfine spin state confined in one-dimensional space. 
In all-optical dipole traps, the three spin substates \F — 
1 , in p = 1), |F = l,m F = 0), and |F = l,m F = -1), 
where mp is the magnetic quantum number, are essen- 
tially free and the spinor nature of the alkali atoms can 
be manifested. In this paper, we discuss the soliton prop- 
erties of this system in the framework of the integrable 
model, which brings us to understand the intrinsic as- 
pects of this complex nonlinear system and provides a 
rigid reference for further numerical and experimental 
investigations. 

The paper is organized as follows. In § 2, the model of 
spinor condensate is presented, and its embedding into 
the integrable equation is reviewed. In § 3, we briefly 
summarize the inverse scattering method, the way to the 
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exact A-soliton solutions of this model, and provide some 
conserved quantities. Next, in § 4, the properties of the 
one-soliton solution classified by the spin states are in- 
vestigated in detail. In particular, we introduce a singlet 
pair density to resolve energetically degenerate solitons 
in the polar state. Then, in § 5, collisional effects between 
solitons in the same or different spin state(s) are exam- 
ined by analyzing the asymptotic forms of their initial 
and final state. As one of the most significant phenomena 
of the two-soliton collisions, spin switching is predicted. 
Section 6 concludes the paper. In Appendix, we give the 
explicit formula of general two-soliton solutions in this 
spinor model. 

2. F — 1 Spinor Condensate 

Atoms in the F = 1 state arc characterized by a 
vectorial field operator with the components subject to 
the hyperfine spin manifold. The three-component field 
* = * , *-i} T , with the superscript T denotes the 
transpose, satisfies the bosonic commutation relations: 

[* a {x, t),ty(x', t)] = 5 a , 5(x - x'). (1) 

Here we assume that the system is one dimensional (ID): 
the trap is suitably anisotropic such that the transverse 
spatial degrees of freedom is factorized from the longi- 
tudinal and all the hyperfine states are in transverse 
ground state. Those quasi-one dimensional regime can 
be achieved experimentally (see the discussion in refs. 
8). The interaction between atoms in the F — 1 hyper- 
fine state is given by 15 ' 16 



V{x 1 - x 2 ) = S(xi - x 2 )(c + c 2 F x ■ F 2 ), 



(2) 



where Fi are the spin operators of two atoms. In this 
expression, 

c = 47rft 2 (ao + 2a 2 )/3m, 

c 2 = 4irh 2 (a 2 - a )/3m, (3) 

where a / are the s-wave scattering lengths for the channel 
of total hyperfine spin / and m is the mass of the atom. 
As a result of the symmetry required for bosonic atoms, 
it can be shown that only even-/ states contribute to 
V(xi — x 2 ). Therefore, in the case of F — 1 atoms, there 
are two scattering channels: / = 2 for the parallel spin 
collision and / = for the antiparallel. According to ref. 
24, the effective ID couplings co, 2 are represented by 

c = c /2a 2 ± , c 2 = c 2 /2a 2 ± , (4) 

where a± is the size of the transverse ground state. 
Thus, the second-quantized Hamiltonian is 

h = J Ax (J^dM ■ d x i> a + c ^¥J>l^ a ^ a 

+ f*i*L/J/3-/a'/J'*/J'*/j), (5) 

where repeated subscripts {a, (3, a', (3' — 1, 0,-1} should 
be summed up. The explicit form of the above Hamilto- 



nian is shown to be 



+ ^i^l!*-!*-! + 2* *o(*I*l + *Li*-l) 



,T,t 



+ (co - c 2 )* t 1 * t _ 1 *_ 1 * 1 + ^* * *o*o 



£o&t,f,ti 



+ C 2 



44*^*0*0 + ^o^O*-!*! 



t,f,t, 



(6) 



where the following expression of spin-1 matrices / 
{f x , P, f z } T is employed, 

, /010\ . / -1 

r = -jR i o i ),f v = As 1 - 1 

VM 1 / ^ 2 \01 



r 




(7) 



In the mean-field theory of BECs, the three-component 
condensate wave function &(x, t) is obtained by 

&(x,t) = (4f(x,t)) 

= {$i(z,*),$o(z,i),$-i(M)} T , (8) 
which is normalized to the total number of atoms At: 



(9) 



The time-evolution of spinor condensate wave function 
<&(x,t) can be derived from the variational principle: 

SEqp 



ihd t $ a (x,t) 



(10) 



where the Gross-Pitaevskii energy functional is given by 

co + c 2 



$* . g $ + 



|$l| 4 +|$- 



+ 2|$ | 2 (|$i| 2 + l^il 2 ) + (co ~ c 2 )|*i| 2 |^-i 



+ f |*ol 4 + c^bVbUn + *S (ii) 

Substituting eq.(ll) into eq.(10), we get a set of equa- 
tions for the spinor condensate wave functions: 

h 2 

m$l = -7^a 2 $! + (Co + C 2 ){|$!| 2 + |$o| 2 }$l 

2m 

+ (Go - c 2 )|*_i| 2 *i + c 2 ^*_^l, 

ifia t $ - -^d 2 $o + (co + c 2 ){|$i| 2 + |*-i| 2 }$ 
+ c |$o| 2 4'o + 2c 2 $*$ 1 $_i, 

h 2 



2m 



9 2 $_i + (co + c 2 ){|$_!i 2 + i* r}*-i 



+ (co - c 2 )|$i| 2 $_i + c 2 <f>l%. 



(12) 



To analyze the dynamical properties of these coupled 
system, we employ the integrable model found in ref. 23. 
We consider the system with the integrable condition 
of coupling constants c = c 2 = — c < 0, cquivalcntly 
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scattering lengths 2a = — a 2 > 0. This situation cor- 
responds to attractive mean-field interaction and ferro- 
magnetic spin-exchange interaction. The effective inter- 
actions between atoms in a BEC can be tuned with the 
so-called Feshbach resonance. 25,26 Moreover, alternative 
techniques, which do not affect the rotational symme- 
try of the internal spin states, such as optically induced 
Feshbach resonance 27-29 are available. 

In the dimensionless form: <fr — > {fa, V%4>o, ^-i} T i 
where time and length are measured respectively in units 
of i — ha±/c and x — hy / a±/2mc 7 we rewrite eqs.(12) 
as follows, 

idtfa = -d 2 ^ - 2{|^| 2 + 2|</, o | 2 }0 1 - 24,*_ 1 4>l 

id t <f> a = ~d 2 x ^ - 2{|0_!| 2 + |^o| 2 + |^| 2 }0 O - 205M_ 1) 



In the following, we shall exploit the exact solution 
technique which enables us to solve the initial value prob- 
lem of the coupled system (14) with (15) and obtain the 
exact soliton solutions. 

3. Inverse Scattering Method 

In this section, some useful results obtained by the 
inverse scattering method (ISM) are briefly summarized. 
We derive an explicit formula for the soliton solution of 
the 2x2 matrix version of NLS equation (14) with eq.(15) 
by considering a reduction of a formula in ref. 30. 

Under the vanishing boundary conditions, we apply 
the ISM to the nonlinear time evolution equation (14) 
associated with the generalized Zakharov-Shabat eigen- 
value problem: 



id t (f>- 



- 2{|<M 2 + 2|^ o | 2 }0_ 1 - 2^. (13) 



These coupled equations are equivalent to a 2 x 2 matrix 
version of nonlinear Schrodinger (NLS) equation: 

id t Q + d 2 x Q + 2QQ t Q = 0, (14) 

with an identification, 



Q = 



I>1 <Po 



(15) 



Since the matrix NLS equation (14) is completely inte- 
grable, 30 the integrability of the reduced equations (13) 
are proved automatically. 23 Remark that the general 
M x L matrix NLS equation is also integrable. 30 It is 
worthy to search other integrable models for higher spin 
case. We will discuss such possibilities in a separate pa- 
per. 

We also mention another reduction: 



<7i 92 




(16) 



This gives rise to intensity-coupled NLS equations known 
as the Manakov model: 31 

i&gi+^«i + 2(|gi| 2 + |<fc| 2 ) gi =0, 

id t q2 + d 2 x q 2 + 2(\ gi \ 2 + \q 2 \ 2 )q 2 = 0, (17) 

which are used to describe the interaction among the 
modes in nonlinear optics, e.g., in the case of bircfrin- 
gent and other two-mode fibers. 32 The extension of the 
Manakov model (17) to the general m-component case 
is straightforward. Recently, ./V-soliton collisions in this 
model have been analyzed in detail. 33 It should be noted 
that by rescaling these two fields q\ , q 2 appropriately, one 
can introduce the coupling coefficients in the nonlinear 
terms. However, possible choices of the coefficients that 
keep the integrability of the system are highly restricted 
according to the Painlcvc analysis. 34 For example, one 
can apply the integrable equations (17) to a binary BEC 9 
only if both the inter-species coupling and intra-species 
coupling coincide. Therefore, in the Manakov model (17), 
two (or more in the general case) components are cou- 
pled equally with each other. On the other hand, in our 
spinor model (13), it is obvious that component and 
±1 components play different roles; this difference will 
lead to new types of the soliton solutions and the unique 
spin dynamics. 



(18) 



1 / k*I 2Q \ *i 

2 ^ -2Q -k*I J { *„ 

Here "Ji and \I/n take their values in 2 x 2 matrices. The 
complex number k is the spectral parameter. / is the 
2x2 unit matrix. The 2x2 matrix Q plays a role as a 
potential function in this linear system. 

The general ./V-soliton solution of eq.(14) with eq.(15) 
is expressed as 30 

/ n ie xi \ 

n 2 e X2 

Q(x,t) = {II^DS- 1 . , (19) 

JV 



\ n N e XN 

where the 2N x 2N matrix S is given by 



N 

E 



^{h + k*)(k 3 + k*) 



(20) 



K) 

1 < i,j < N. 

Here we have introduced the following parametcrizations: 



n, = 



Xj = Xj{x, t) = kjX + \k)t - ej. (21) 

We explain notations and their significance. The 2x2 
matrices H,- normalized to unity in a sense of the square 



norm, 



in 



jl|2 



2\aj\ 2 



i, 



(22) 



must take the same form as Q from their definition. We 
call them "polarization matrices" as usual in the Man- 
akov model. In the spinor model, the polarization matri- 
ces determine both the populations of three components 
{1, 0, —1} within each soliton and the relative phases be- 
tween them. The complex constants kj denote discrete 
eigenvalues, each of which determines a bound state by 
the potential Q. ej are real constants which can be used 
to tune the initial displacements of solitons. It is worth 
noting that all x and t dependence is only through the 
variables \j{x,t). As we shall see in section 4, the real 
part of \j(x,t) represents the coordinate for observing 
soliton- j's envelope while the imaginary part of it rep- 
resents the coordinate for observing soliton-j's carrier 
waves. 
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Equation (14) is a completely integrable system whose 
initial value problems can be solved via, for example, 
the ISM. The existence of the r-matrix for this system 
guarantees the existence of an infinite number of conser- 
vation laws which restrict the dynamics of the system in 
an essential way. Here we show explicit forms of some 
conserved quantities, i.e., total number, total spin, total 



momentum and total energy, 

total number: Nt = J dxn(x,t); (23) 

n(x, i) = *t . $ = tr{Q^Q}. (24) 

total spin: F T = J dxf(x,t); (25) 

f(x,t) = * t -/-* = tr{Q t crQ}. (26) 

total momentum: P T = J dxp(x,t); (27) 

p(x, t) = -ift*t . dx <£ = _yj . tr{QtQ x }. (28) 

total energy: Et = J dxe(x, t) ; (29) 

e(x, t) = ^-d x ^ ■ d x & - | [n(x, t) 2 + f(x, t) 2 ] 

= C -tr{QtQ x _QtQgtQ}. (3 ) 



Here tr{-} denotes the matrix trace. In eq.(26), the ex- 
pression (7) and the Pauli matrices <x = (<T x ,<r y , a z ) T are 
used. These conservation laws can be proved by a direct 
calculation. For a complete set of the conserved densities 
and the recursion formula for them, see discussion in ref. 
30. 

4. Spin State of One-Soliton Solution 

In this section, we discuss one-soliton solution and the 
classification of its spin states in the spinor model. If we 
set N — 1 in the formula (19), we obtain the one-soliton 
solution: 

_ e* / (i + 7 *e 2xR+p dctn a - a*e 2xR+ PdetU \ 
Q ~ dctS \ a- a*e 2xR+ PdetU 7 + /3*e 2xR+p dctn ) ' 

(31) 

where 

dctS* = 1 + c 2 ^ +p + e 4xR - +2 "|dctn| 2 , (32) 

e p / 2 = — — , U=( P °), (33) 

XR = Xr(x, t) = k R (x - 2kit) - e, (34) 
Xi = Xi(x, t) = hx + (/c 2 - k 2 )t. (35) 

We have omitted the subscripts of the soliton number. 
Here and henceforth, the subscripts R and I denote real 
and imaginary parts, respectively. Throughout this sec- 
tion, we set k R > without loss of generality. In a dif- 
ferent form: 

Il e -(xn+p/2) + (o-s/nto-s/) e^+^deffl ; 
= 2fcR e -( 2 XH+p) + l + e 2 XH+P|detn| 2 ^ ' (36) 



we can make out the significance of each parame- 
ter/coordinate as follows, 

n : polarization matrix of soliton 

: amplitude of soliton 

2ki : velocity of soliton's envelope 

Xr : coordinate for observing soliton's envelope 

Xi : coordinate for observing soliton's carrier waves. 

A few comments will be made on those points pertinent 
to this list. We use the term "amplitude" to indicate the 
peak(s) height of soliton's envelope. Actual amplitude 
should be represented as k R multiplied by a factor from 
1 to \/2 which is determined by the type of polarization 
matrices. The explicit form will be shown later. As men- 
tioned before, soliton's motion depends on both x and 
t via variables xr and Xii from which we can see the 
meaning of velocity of soliton. 

From a total spin conservation, one-soliton solution 
can be classified by the spin states. We shall show 
that the only two spin states are allowable, i.e., Ft = 
(0, 0, 0) T for deffl ^ 0, and |F T | = N T for deffl = 0. 

Substituting eqs.(31)-(35) into eq.(26), we obtain the 
local spin density of the one-soliton solution: 

p 2XR 

f(M) = (deW (1 ~ c4XR+2p l dcffl l 2 ) ^nVn}. (37) 

We also give the explicit form of the number density: 

n{x,t) = {1 + (4c 2 »+" + e 4xR+2 ") |detn| 2 } . 

(38) 

To clarify the physical meaning of det n, we define here 
another important local density as 

G(x,t) = ^ - 2$i$_i = -2dctQ. (39) 

This quantity measures the formation of singlet pairs. 
Note that these "pairs" are distinguished from Cooper 
pairs of electrons or those of 3 He owing to the different 
statistical properties of ingredient particles. 18 ' 20 Since 
0(x,t) does not contribute to the magnetization of the 
soliton, it is invariant under any spin rotation. 19 As far as 
ground state properties are concerned, it is not necessary 
to introduce <d(x,t) for a system of spin-1 bosons, while 
a counterpart to eq.(39) plays a crucial role for spin-2 
case. 18-20 However, as we shall show later, it is useful to 
characterize solitons within energy degenerated states. 

In the case of the one-soliton solution (31), the singlet 
pair density is proportional to the determinant of the 
polarization matrix n, 

e 2x 

G(x, t) = -2- — - det n. (40) 
det b 

This suggests that det n represents the magnitude of the 
singlet pairs locally. For the general case, however, this 
singlet pair density can vary after each collision of soli- 
tons and is not the conserved density. The detail will be 
discussed at the end of this section and also in §5.2. Now, 
we classify spin states of the one-soliton solution into two 
distinct cases depending on the values of det n. 
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4-1 Ferromagnetic state 

Let detLT = 0, then eq.(31) becomes a simple form: 



Q = fc R sech( X R + p/2)ne i; 



(41) 



Now all of m f = 0, ±1 components share the same wave 
function. Their distribution in the internal state reflects 
the elements of the polarization matrix II directly. One 
can clearly see the meaning of each parameter listed 
above. By definition, the singlet pair density (40) van- 
ishes everywhere. Thus, this type of soliton belongs to 
the ferromagnetic state and will be referred to as a ferro- 
magnetic soliton. The total number of atoms is obtained 
by integrating eq.(38) as Nt — 2fc R . The total spin (25) 
becomes 



F T = 2fc F 



2Re{a*(/3 + 7 )} 
-2Im{a*(/3- 7 )} 
l/3| 2 -M 2 



|F T I = Nrr 



(42) 



which is connected to F T = 2fc R (0,0, 1) T through a 
gauge transformation and a spin rotation. 

Next, we calculate the total momentum and the total 
energy of the ferromagnetic soliton. Substituting cq.(41) 
into eqs.(27), (29), and using detll = 0, we obtain 



P f 



-ihk 



R J dx {ifci - k R tanh(xR + p/2)} 
x sech 2 (xR + p/2) 



N T hh, 



(43) 



e t = ck R J Ax [ k R {2 tanh 2 (xR + p/2) - 1} + kf 



x sech 2 (xR + p/2) 

2 



(44) 



respectively. 



4-2 Polar state 

If detll 7^ 0, the local spin density has one node, i.e., 
f(xo,t) = at a point: 



1 / 4fc 2 



2e 



(45) 



for each moment t. Setting x' = x — x and A 1 
2 1 dctil|, wc get 

2Rc{a*(/3 + 7 )} 



4fc R ylsinh(2fc R x') 
[cosh(2fc R x') + A] 



-2Im{a*(/3- 7 )} 

lfl|2 u,|2 



(46) 

Since each component of the local spin density is an odd 
function of x', its average value is zero, 



dx' f (x') = ( 



(47) 



This implies that this type of soliton, on the average, be- 
longs to the polar state. 15 Let us also rewrite the number 
density (38) as 

4fc R [Acosh2fc R x' + 1] 



n(x') 



[cosh2k R x' + AY 



(48) 



To elaborate on this type of soliton, we further divide 
into two cases. 

(i) A- 1 = 2| dctn| = 1 {af3* + a* 7 =0). 

Under this constraint, we find the local spin (46) itself 
vanishes everywhere. Solitons in this state possess the 
symmetry of polar state locally. We, therefore, refer to 
only those solitons as polar solitons in this paper. Con- 
sidering eq.(36) with the above condition, we recover a 
normal sech-type soliton solution: 

V / 2fc R sech(/c R a;')ne ixi . 



Q 



(49) 



Note that the amplitude of soliton is different from that 
of the ferromagnetic soliton, which leads to a relation 
between the total number and the spectral parameter as 
Nt = 4fc R . The total momentum and the total energy 
are given by P| = N T hk h E T = N T c (fcj 2 - fc|/3), re- 
spectively. The difference between ferromagnetic soliton 
energy and polar soliton energy with the same number 
of atoms Nt is 



16 



<0, 



(50) 



which is a natural consequence of the ferromagnetic 
interaction, i.e., c 2 < 0. 

(ii) A- 1 = 2|detn| < 1. 

In this case, the local spin retains nonzero value, al- 
though the average spin amounts to be zero. The den- 
sity profile has a complicated form as eq.(48). Note that 
when A > 2, a peak of the density splits into two (Fig.l) 
due to different density profiles of mp = 0, ±1 compo- 
nents. For a large value of A, viz., when detLT gets close 
to zero, such twin peaks separate away. In consequence, 
they behave as if a pair of two distinct ferromagnetic soli- 
tons with antiparallel spins, traveling in parallel with the 
same velocity and the amplitudes half as much as that 
of the polar soliton (A = 1) in the density profile (see 
the inset of Fig.l and Fig. 2). Hence, solitons of this type 
will be referred to as split solitons. The total number is 
the same as (i) case, A^t = 4fc R , which can be shown 
by changing a variable y = tanh(fc R x') and evaluating a 




-6 -4 -2 



Fig. 1. The density profiles of eq. (48). We set Air = 0.5, and 
A = 1 (solid line), 2 (dashed line), 5 (dash-dot line), 20 (dotted 
line). The inset shows a split soliton for A = 10 4 , consisting of 
two ferromagnetic like solitons with the same velocity. 
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parameters to regions: 





5 10 15 * 

Fig. 2. The density profiles of eq.(48) (solid line) for = 0.5 
and A = 10 4 , and the three components, mp = (dashed line), 
mp = 1 (dotted line) and mp = —I (dash-dot line) are shown 
simultaneously. 



definite integral: 



i , (A+l) + y 2 (A-l) 
i [{A+l)-y\A-l)f 



1. 



(51) 



By use of the similar definite integrals, the total mo- 
mentum and the total energy are shown to be the same 
values as those in the case (i): Prf, = N T hki, E^ — 



N T c (Af 



fcp/3) . This degeneracy is ascribed to the inte- 



grable condition for the coupling constants, i.e., Co = c 2 . 

Comparing case (i) with case (ii) , we find that a variety 
of dissimilar shaped solitons are degenerated in the polar 
state. To characterize them, we can use, instead of A 
itself, a physical quantity defined as 



S = / dx\e(x, t)\ 



2 tan 



= N'T 



1 (V*5^ 



1 



(52) 



which is a monotone decreasing function of A e [1, oo); 
the maximum value, Nt, at A = 1 (polar soliton) and 
limiting to in A — > oo (ferromagnetic soliton). In this 
sense, S has the meaning as the "total singlet pairs" of 
the whole system. As noted above, S is not the conserved 
quantity in general (N > 2); all the conserved densities 
should be expressed by the matrix trace of products of 
Qt, Q and their derivatives 30 as eqs.(24), (26), (28), and 
(30) while \Q{x, t)\ is not. Nevertheless, S can be used 
to label solitons in the polar state because it dose not 
change in the meanwhile prior to the subsequent colli- 
sion. 

5. Two-Soliton Collision and Scattering Law 

In this section, we analyze two-soliton collisions in the 
spinor model. The two-soliton solution can be obtained 
by setting N = 2 in eq.(19). The derivation is straight- 
forward but rather lengthy. We give a general formula 
of the two-soliton solution in Appendix and, here, com- 
pute asymptotic forms of specific two-soliton solutions as 
t — ► =F°o, which define the collision laws of two-soliton in 
the spinor model. For simplicity, we restrict the spectral 



fciR > 0, fc 2R < 0, 
fcn < 0, k 2 i > 0. 



(53) 



Under the conditions, we calculate the asymptotic forms 
in the final state (t — ► oo) from those in the initial 
state (t — ► — oo). Since each soliton's envelope is located 
around x ~ 2kjjt, soliton- 1 and soliton-2 are initially 
isolated at x — > ±00, and then, travel to the opposite di- 
rection at a velocity of 2k\\ and 2fc 2 i, respectively. After 
a head-on collision, they pass through without changing 
their velocities and arrive at x — > =Foo in the final state. 
Collisional effects appear not only as usual phase shifts 
of solitons but also as a rotation of their polarization. 

According to the classification of one-soliton solu- 
tions in the previous section, we choose the following 
three cases: (1) polar-polar solitons collision, (2) polar- 
ferromagnetic solitons collision, and (3) ferromagnetic- 
ferromagnetic solitons collision. As we shall see later, the 
polar soliton dose not affect the polarization of the other 
solitons apart from the total phase factor. On the other 
hand, ferromagnetic solitons can 'rotate' their partners' 
polarization, which allows for switching among the inter- 
nal states. 

5.1 Polar-polar solitons collision 

We first deal with a collision between two polar soli- 
tons defined by kj and II, (j = 1, 2) with the conditions 
(53) and |detlli| = |detn 2 | = 1/2. In the asymptotic 
regions, we can consider each soliton separately. Thus, 
the initial state is given by the sum of two polar solitons 
as 

Q * Q'T + Q 2 n , (54) 

where the asymptotic form of soliton- j (j = 1, 2) is 

Qf = \/2fc jR sech( Xj R + />, - > ;ll.,<- ix . (55) 

These can be proved by taking the limit %2R — * —00 
with keeping \m finite and, vice versa, Xir - * with 
X2R fixed. Phase factors which come from the values of 
I dct IL, I arc absorbed by the arbitrary constants ej inside 
XjR- In the final state, the opposite limit X2R — * 00 with 
keeping xir finite and xir ~ > 00 with |x2r| < 00 yield 



Q - Q? n 



where 
with 



V2fc iR sech (xjr + Pj/2 



2 In 



ki 



o~i = 2arg 



fci + k* 
k\ — . 



fci + k% 



0-2 = 2 arg 



k 2 - fci 



(56) 



(57) 



(58) 



•(59) 



Equations (55) and (57) are the same form as polar 
one-soliton solution (49). Collisional effects appear only 
in the position shift (58) and the phase shifts (59). Thus, 
the partial number Nj, spin Fj, momentum Pj, and en- 
ergy Ej are defined for the asymptotic form of soliton- 
j and calculated in the same manner as the previous 
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section. The integrals of motion are represented by the 
sum of those quantities for each soliton. Moreover, we 
can prove Nj — 4|fcj R |, \Fj\ — 0, Pj = Njhkji, and 
Ej = Njc(k^ — kj R /3) themselves are conserved through 
the collision. In this sense, the polar-polar collision is 
basically the same as that of the single-component NLS 
equation. 

5.2 Polar- ferromagnetic solitons collision 

Under the condition (53), we set soliton 1 to be polar 
soliton (|detlli| = 1/2), and soliton 2 to be ferromag- 
netic soliton (| det II2 1 = 0). Then, the initial state is 
represented by eq.(54) with 

Q[ n = V2fc 1R sech( XlR + p 1 /2)n 1 e i X", 

Q 2 n = fe R sech( X2R + p 2 /2)n 2 c^ 21 . (60) 

The final state is given by eq.(56) with 

Qf = 2fc 1R e i * 11 

n lC -(xiR+Pi/2+5) + (^nt^^etn^xiR+Pi/a+.s 

X c -(2 XlR + Pl +2(5) + l + c 2 X ir+Pi+2<5| detflil 2 ' 

QT = fc 2 Rsech( X 2R + P2/2 + r) n 2 e^ M \ (61) 
Here we have defined 

2 



e 25 = 



ki - k 2 



k\ + k 



2 



(riiiT 2 )| 2 }, 

fii = e-^n! - (n in tn 2 + n 2 nt ni ) 



\h + k^lh - k 2 \ 2 

k 2 + k$ 



+ 



(63) 



fci + ft 2 

and also used eqs.(58), (59). Normalization of the new 
polarization matrix (63) turns to be unity, 1 1 XT x 1 1 2 = 1- 
The determinant of it becomes 

2 



det IIi = c 



fci - k 2 



detni. 



(64) 



We can see clearly that the initial polar soliton breaks 
into a split type (A\ = (2| detlli|) _1 > 1) after the colli- 







sion with a ferromagnetic one. Only when |tr(nill|) 
(in this case, the spinor wavcfunctions of two initial soli- 
tons are orthogonal), we have A\ = 1. Then, cqs.(60) are 
reduced to 

Qf - V2fc 1R sech( XlR + p 1 /2 + r)n 1 c i fe I+CTl ), 

Qt n = fe R sech ( X2R + p 2 /2 + r) H 2 e^ M \ (65) 

which means that the polar soliton keeps its shape after 
the collision and shows no mixing among the internal 
states except for the total phase shift. On the other hand, 
because of the total spin conservation, the ferromagnetic 
soliton always retains its polarization matrix and shows 
only the position and phase shifts similar to those of the 
polar-polar case. 

In Fig. 3, we have density plots of a polar- 
ferromagnetic collision with the parameters shown in the 
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Fig. 3. Density plots of |<^o| 2 (top) and |</>±i| 2 (bottom) for a 
polar-ferromagnetic collision. Soliton 1 (left mover) is a polar 
soliton and soliton 2 (right mover) is a ferromagnetic soliton. The 
parameters used here are fci = 0.25 — 0.25i, fc2 = —0.5 + 0.25i, 
ai = 0, f3 1 = 7i = 1/V2, oc 2 = #2 = 72 = 1/2. 
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Fig. 4. Density plots of \4>o\ 2 (dots) and |</>±i| 2 (solid line) in the 
initial state (bottom) and the final state (top) of Fig. 3. 



caption. These pictures correspond to each component 
of the exact two-soliton solution (A- 9) for one collisional 
run. For simplicity, we choose the parameters to have 
|0i I = l^-i I • Figure 4 shows the overlap of |</>o| 2 and 
|0±i| 2 in the initial state (60) and in the final state 
(61). The polar soliton (soliton 1) initially prepared in 
tuf = ±1 are switched into a soliton with a large popu- 
lation in rap = and the remnant of mj? = ±1 after the 
collision. Through the collision, the ferromagnetic soliton 
(soliton 2) plays only a switcher, showing no mixing in 
the internal state of itself outside the collisional region, 
as clearly seen in eq.(61). In general, this kind of a drastic 
internal shift of polar soliton is likely observed for large 
values of |tr(lliil|) | which appears in eqs.(62), (63). Al- 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



though all the conserved quantities such as the number 
of particles and the averaged spin of individual solitons 
are invariant during this type of collision, the fraction of 
each component can vary not only in each soliton level 
but also in the total after the collision. This contrasts to 
an intensity coupled multi-component system in which 
the total distribution among all components is invariant 
throughout soliton collisions. 

5.3 Ferromagnetic-ferromagnetic solitons collision 

As an instructive example, we already discussed the 
collision between two ferromagnetic solitons (det II i = 
det n 2 = 0) in ref. 23. In this subsection, we recapitulate 
it for completeness of this paper and extend the analysis 
of the spin-precession dynamics to the general case. 

The asymptotic forms are obtained for the initial state, 
Q ~ + Q 2 n where 

Qf = fc, R sech( X ,R, + 1>; 11,,^ '. 

and for the final state, Q ~ Q\ n + Q 2 n where 

Qf = fc jR sech ( XjR + p^l + s) flje^ . 

Here we have defined 

(fci+fcrxjfc+fcs). 



\ni 1 



Ifcl 



and, for = (1,2) or (2,1) 

( Uj - kl+k - 
kj + k'l 



|tr(n 1 n+) 



(66) 
(67) 
(68) 



n, = e - 



+ 



h + k* 



tr^-nDn, , 



(69) 



which are shown to be normalized in unity, 1 1 n ^ 1 1 2 = 1- 
Each polarization matrix IX,- of a ferromagnetic soliton 
can be expressed by three real variables Tj , 9j , ipj as 



n, = e" 



cos z -fe 1Vj cos -f sin -± 

2 2 2 

cos — sin — sin — e Vj 

2 2 2 



(70) 



In this expression, the polarization matrices in the initial 
state II, and in the final state IL, are given by 

(71) 
(72) 



n, = e^ Uj 



n, = e" 



S+1T, 



Uj ■ u, 



where, with = (1,2), (2,1), 
/ 9, 



u, 



V 



cos -re 1 



sin — e 2 
2 



k, 



(73) 



(74) 



k j + k *i 

This defines the collision law for the ferromagnetic- 
ferromagnetic soliton collision. Similar collision laws for 
the matrix KdV equation can be found in ref. 35. 



We calculate the spins for each soliton to investigate 
their collision. It will be interpreted as a spin preces- 
sion around the total spin. In the initial state, following 
eq.(42), we have the spin of soliton- j as 



F, = 2\k. 



sin 9j cos ifj 
K j sin 9j sin ipj 
cos 9j 



(75) 



Thanks to the scattering law (74), the final state spins 
can be obtained through Fi ; 2 by 



1 _ 2fc; R (fci R + k 2 R, . F 



\ki+k*V ' '" ' 



+ \hTW {F]XFl) 

+ | eS _ 1 + 2fc jR (fci R + h R ) \ 



\h+ktf J 



where 



e s = 1 + 



|Fi|-|F 2 |) 2 



(76) 



(77) 



4|fc!+fc 2 *| 2 

The conserved total spin, F T = Fi + F 2 = Fi + F 2 , is 
given by 

(2 1 fci R I sin 9\ cos tpi + 2| fc 2R | sin 62 cos ip 2 
2|fci R | sin 9\ sin<pi + 2|fc 2R | sin# 2 sin<p 2 
2|fci R | cos 6»i + 2|fc 2R | cos 2 

(78) 

Considering spin rotation around the total spin Ft, 
we can find 'rotated spin' as 

Ff - r 1 (i P )h- 1 (9)f( l j)h(9)f( l p)F j , (79) 

where 

cos cp sin ip 
— sin ip cos ip 
1 



m 

h(9) 



cos 9 — sin 9 

1 
sin 9 cos 9 



(80) 



and 

COSt/3 



= , sin 93 



y/(F§)* + (F|) 2 



cos 9 



sm( 



V(^) 2 + (^|) 2 ' 
(81) 



Then, 



-17 rot ^ 

F ^ =2 



{(1 + Aji) + (1- Aj-Ocosw)}^ 
(F, x Fj) 

+ {(1 + A ji ) - (l + A j j)cosw}Fi 



,(82) 



where 
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CO/ji 




Fig. 5. Velocity dependence of the rotational angle in spin pre- 
cession for the different initial relative angles between two spins, 
T = 1 (solid line), 0.5 (dashed line), 0.0157-7T (dash-dot line) and 
(dotted line). 



The rotation angle uj is determined by settin; 
through eqs. (76) and (82), 



-prot 



COSW = 



sinw 



4(fcn - fc 2 i) 2 



4(fcn - M 2 + |Ft| 2: 
4(fc 2 i-fcn)|F T | 



(84) 



(85) 



4(fcn - fc 2 i) 2 + |F T | 2 " 

For the case that the amplitude and velocity of two 
ferromagnetic solitons are same, | | = |/c 2R | = N T /A 
and | fen | = |fc 2 i| = fci, the final state spins Fj are given 

by 

(Fj x F,) 



F, = cos- ^F, - ! - sin - ]~ r ' ; +sin 2 



(86) 



where (J, I) — (1, 2), (2, 1). The rotation angle uj depends 
only on the ratio ki / kn and the magnitude of the nor- 
malized total spin T = \Ft/Nt\ as 

-l/2\ 



uj = 2 arccos 



1 + 



ki 



(87) 



The principal value should be taken for the arccosine 
function: < arccos x < ir. 

Setting ki ^> k# in eq. (87), one gets the small rota- 
tion angle, uj ~ 0. In the opposite case, k\ <C fc R , each 
spin of two colliding solitons almost reverses its orienta- 
tion, uj ~ 7r. Recall that k\ is the speed of soliton. We 
can understand these phenomena since a slower soliton 
spends the longer time inside the collisional region. Fig- 
ure 5 shows the velocity dependence of the rotation angle 
for various initial normalized spins. When T = 0, which 
corresponds to the case of antiparallel spin collision, the 
spin precession can not occur as shown by the dotted line 
in Fig. 5. 

In Fig. 6-Fig. 8, we give examples of this type of col- 
lisions changing k\, with the other conditions fixed, to 
illustrate the velocity dependence. The initial normal- 
ized spin for the parameter set given in the captions is 
T = 0.5. The rotation angles are uj ~ 0.2n, 0.5tt and 
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Fig. 6. Density plots of (a) \(po\ 2 , (b) |</>i| 2 and (c) |0-i| 2 for a 
fast ferromagnetic-ferromagnetic collision. The parameters used 
here arc fci = 0.5 - 0.75i, k 2 = -0.5 + 0.75i, «i = 4/17, ft = 
16/17, 7i = 1/17, a 2 = 4/17, #2 = 1/17, 72 = 16/17. 
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Fig. 7. Density plots of (a) |</> | 2 , (b) Iftl 2 and (c) \4>-i\ 2 for 
a medium speed ferromagnetic-ferromagnetic collision. The pa- 
rameters are the same as those of Fig. 6 except for kn = —0.25, 
fc 2 i = 0.25. 
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Fig. 8. Density plots of (a) |0 O | 2 , (b) |</>i| 2 and (c) |0-i| 2 for a 
slow ferromagnetic-ferromagnetic collision. The parameters are 
the same as those of Fig. 6 except for fen = —0.05, k 2 i = 0.05. 



0.97T for Fig. 6, Fig. 7 and Fig. 8, respectively. The in- 
ternal shift <fii — ► 4>-i, and vice versa, gradually increase 
by slowing down the velocity of the solitons. 
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6. Conclusions and Discussions 

In the present paper, we have studied the soliton prop- 
erties in F = 1 spinor Bose-Einstein condensates on the 
basis of the integrable model introduced in ref. 23. 

The complete classification of the one-soliton solution 
has been carried out. There exist two distinct spin states, 
ferromagnetic, |Ft| = Nt and polar, |Ft| = 0. In the fer- 
romagnetic state, the spatial part and the spinor part of 
the soliton decouple (ferromagnetic soliton). In the po- 
lar state, dissimilar shaped solitons which we call polar 
soliton for f(x) = and split soliton otherwise are ener- 
getically degenerate. The polar soliton has one peak and 
the space-spinor decoupling holds. On the other hand, a 
split soliton consists of twin peaks and the profiles of the 
three components are different. Changing the polariza- 
tion parameters, one may control the distance between 
these peaks continuously. 

These properties remind us of the Sasa-Satsuma 
higher-order nonlinear Schrodinger (HNLS) equation, 
which is a unique soliton equation generating solitons 
with two peaks. In fact, the spinor model and the Sasa- 
Satsuma model, though the former is three-component 
and the latter is one-component, have exactly the same 
form of the envelope, i.e., the density (48). (See also 
eq.(52a) in ref. 36.) The relation between the multi- 
component NLS equation and the one-component HNLS 
equation is interesting and remains as a future problem. 

We have also shown explicit two-soliton solutions 
which rule collisional phenomena of the multiple solitons. 
Specifying the initial conditions, we have demonstrated 
two-soliton collisions in three characteristic cases: polar- 
polar, polar-ferromagnetic, ferromagnetic-ferromagnetic. 
In their collisions, the polar soliton is always "passive" 
which means that it can not rotate its partner's polar- 
ization while the ferromagnetic soliton does. Thus, in the 
polar-ferromagnetic collision, one can use the polar soli- 
ton as a signal and ferromagnetic soliton as a switch to 
realize a coherent matter- wave switching device. Colli- 
sion of two ferromagnetic solitons can be interpreted as 
the spin precession around the total spin. The rotation 
angle depends on the total spin, amplitude and velocity 
of the solitons. Only varying the velocity induces drastic 
change of the population shifts among the components. 

Remark that in current experimental setups, 22 the ini- 
tial population among the three components can be con- 
trolled by magnetic field gradients during preparing con- 
densates. The components overlapping in solitons can be 
spatially separated by applying a weak Stern-Gerlach 
field before a time-of-flight measurement. We hope that 
those phenomena predicted in this work are observed in 
experiments and open up a variety of applications in co- 
herent atom transport, and quantum information. 
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Appendix: Two-Soliton Solution of Spinor Model 

In this appendix, we write out general two-soliton solution given by (19) with N — 2. According to (20), the 
matrix S in this case takes the form: 



+ 



M- 



M< 



ki + kX 



12 c xi+x? 



M 



h + fc* 



ki + k. 



k 2 + k{ 



11 e xi+xl _|_ 



M 



k\ 

H2_ e X2+xt + M22 c y2+y; j + _^_ e X2+xl 



Ah 



k 2 + k{ 



k 2 + kl 
M 22 

k 2 + k*' 



I 2 . gXl+X2 



aX2+X 2 



(A-l) 



where 2x2 matrices Mjj are 



Mi. 



Kij 



(A-2) 



with 



a t a* + faff 
ki ~t~ k; 



Ay — 



a; 7* + 

ki foj 



7ia* + ai/3* 

&i H~ fcj 



Gnu* + 7,7* 
&i H~ 



(A-3) 



Here the elements of the polarization matrices aij, f3j, jj, and the spectral parameters kj are determined by the 
initial state (t — ► — oo). Recall that all x and i dependence is only through the complex coordinates Xj = Xj 
The determinant of (A-l) is 



detS = l + J2 A 3i eX3+x ' + J2 B ljkl e x ' +X ^ +Xk+X ' 



3,1 



i,j,k,l 



C l0k ie 2Xt+2x ^ +Xk+x ' +L» c 2 (xi+xJ+X2+x2) ; 

i,j,k,l 



(A4) 



where 



tiMji 



(A-5) 



Bijki 



ijkl 



D 



det 



Kij Xij 



+ \ {tr(M y M fcJ ) - tr(M«My)} 



(fc 4 + fc*)(fc fe + fc ; *) 
1 



fc 4 + fc* (*< + fc*) (fe fc + fcf ) {h + k*){k k + k*) 



x < det 



Mij 




f ^_ 

~ \ (fci + fcf)(fc 2 



Then, the two-soliton solution is given by 



det 



Afn 
M 2i 



M12 

M 22 



My 



(A-6) 



(A-7) 
(A8) 



Q 



det 5 
1 

det 5 



ii + ,52i)nie Xl + (Sia + S 22 )n 2 e> 



{Sll + S£)/3i + {S£ + SH) ai {S\\ + S£) ai + (S$ + S 2 12 ) 7 i 
(SR + S*l)h + (£i 2 i + 4 2 i 2 )«i + Sl\) ai + (Sfi + 4 2 i 2 )7i 



+ 



(Si?, + S%)/32 + (SH + S%)a 2 (S£ + S£)a 2 + (SH + 5 2 12 ) 72 
(S% + SM)P 2 + {Sit + SE)a 2 {SH + SH)a 2 + {S™ + S 2 22 ) 72 



(A9) 



Here we use the tilde to denote cofactors. The cofactor SfJ is obtained by deleting the (2i + k — 2)-th row and the 
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(2j + I — 2)-th column in the determinant of S and multiplying it by (— l) k+l . For instance, the cofactors S\ l j read 

vij . trM 2 



Si 1 ? 



j=l,2 

+ E 



3 xi+x,- 



fci + fc* 



fc 2 + k 



(fc 2 + fc*) 2 



«2j^lj - A2jMl3 r Xl+X2+2 X : + «21^22 - A 2 1^22 + «22^21 - A 22 ^21 r2y2+x * +X% 



3 tt a {k 1 + k*){k 2 + k*) (fc 2 + fc*)(fc 2 + fc*) 

«22^11 - A 2 i,Ui 2 + ^n^ 22 - ^12^21 , «21^12 - A22MH + ^12^21 - ^11^22 



+ 



(fci + fc*)(fc 2 + fc|) 
|fci-fc 2 | 2 



*|2 E 



^ 2 = E 



J'=l,2 
+ 



(fcj + fcf)(As2 + + fc|| 2 ^ 2 fc 2 + fc; 

. ( ,M < .... — :.„. _- 



(fcl+fc 2 *)(fc 2 +fc*) 

^11 Mij ^12 
det I A 2 i n 2 j A 22 

^21 M2j ^22 



e Xl+X2+xI+X2 



e Xl+2x2+xI+X2+Xj (A-10) 



U r xi+y_: , K ijA2j ^gjAij y 1+ y 2+2y ; 

fci + fc; (h + k*)(k 2 + k*) 



{K12A21 — K22A11 + Al2^21 — An^22 KI1A22 — K2lAl2 + An^22 — Ai2^21 1 Xl+X2+X * +X * 
(fci+fc*)(fc 2 + fc 2 *) (h + k*)(k 2 + kl) 



\ki 



k 2 \ 2 



^|2 E 



(fci + fc*)(fc 2 + fc^lfci + fc 2 *l 2 ^ 2 k 2 + k* 



All Klj A12 

det I A 2 i K2j A 22 

^21 M2j ^22 



7' 



e Xl+2X2+xI+X2+Xj (ATI) 



511 

D 12 



E 

j=l,2 



K U e Xl+X* + 



det M 



k 2 + k 



(h + k*)(k 2 + k*) 



3 2 X i+2 X * 



E 



Klj^j ~ Aij// 2 j cyi +y 2 +2 y ; _ ^11^22 ~ AllMgj + «12^21 - A 12 ^21 ^ t +y 2 + xr +X^ 

f 2 (fe + fc;)(fc2 + A;) (/c 2 + fci*)(fc 2 + fc 2 *) 



K\ 2 V\\ — A11/I12 K\\V\2 — A12MII 



(fci + fc*)(fc 2 + fc 2 *) (fci + fc2)(fc 2 +A;*) 



,2xi+xJ+x 2 



|*1 



E 



(fci + fc?)(fc 2 + + fc 2 T ^ 2 k 2 + k* 



det 



A 



12 



K\j An 

My ^11 ^12 

M2j ^21 ^22 



3 2xi+X2+X^+X2+Xj 



(A-12) 



M2 
'12 



E 

3=1,2 



fc 2 + k* 



e Xl+X* + . Kl J A2 J . . K 2j A lj _ C y 1+X2 +2 X : 



(k 2 +k*)(k 2 +k*) 



\\\v\2 - Aigj^ii Aig^ii - Anz/ia ) 2 X1+X t+ X * 



(fci + • 



(fci + fc2)(fc 2 + 



K11A22 - ^2 2 An + K12A21 - K21A12 Y1+y2+Y r+ y ; 

(fc 2 +fc 1 *)(fc 2 +fc 2 *) 



+ 



|*i - k 2 



(fci + fc*)(fc 2 + k\)\k x + fc. 



*|2 E 



= 1,2 k 2 +k l 



det 





An 


A12 


My 




^12 


«2j 


A21 


A22 



3 2xi +X2 +X* +X*2 +X j 



(A-13) 



The determinant of matrix £ (A-4) is a polynomial in e Xj and e x ^ of degree 8, and the cofactors are degree 6, 
giving rise to a variety of the scattering processes. This property contrasts to the two-soliton solution of the general 
m-component Manakov model, in which the determinant of S and the cofactors are degree 4 and 2, respectively. If 
one set detllj = (j = 1,2), corresponding to the ferromagnetic state, the determinants in the two higher degree's 
coefficients vanish, and the expansions of det S and S^J terminate at degree 4 and 2, respectively. 



